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EFFECT  OF  CHARGE  RADIUS  ON  DETONATION  VELOCITY 


R.  B.  Parlin  and  D.  W.  Robinson 


Abstract 


The  derivation  of  the  detonation  velocity  of  an  ideal  wave  in 
an  imperfect  gas  obeying  the  equation  of  state  of  the  forny 

Vpv  •=  nRTef/^. 

yis  given.  Assuming  the  curved-front  theory,  the  effect  of  radial 
expansion  on  the  velocity  of  propagation  of  a  non-ideal  detonation 
wave  is  treated  using  the  above  equation  of  state.  The  significant 
approximations  involved  are  discussed  and  the  final  equations  are 
given  in  a  form  easily  applied  to  experimental  data.  The  results  of 
this  study  suggest  that  longer  reaction  zone  lengths  than  those 
previously  obtained  using  a  constant  co-volume  may  be  appropriate,  a 


EFFECT  OF  CHARGE  RADIUS  ON  DETONATION  VELOCITY 


R.  B.  Parlin  and  D.  W.  Robinson 

1.  A  qualitative  discussion  of  the  curved-front  theory  is  given 
in  Chem.  Rev.  45,  69  (1949)  (Eyring,  Duffy,  Powell  and  Parlin)  and 

in  Technical  Report  No.  I  (January  1,  1952)  (Parlin,  Thorne  and 
Robinson).  Reference  for  complete  details  of  this  theory  should  be 
made  to  these  reports.  Briefly,  however,  the  theory  explains  the 
observed  decrease  in  detonation  velocity  arising  from  a  reduction 
in  the  radius  of  the  charge.  Consider  a  detonation  wave  propagating 
along  a  charge  of  radius  R  and  of  infinite  length.  Due  to  the  ex¬ 
panding  gas  front,  a  rarefaction  wave  is  sent  into  the  reaction 
zone.  Such  a  rarefaction  wave  will  overtake  the  front  of  the  wave, 
giving  rise  to  a  reduction  in  its  velocity  and  a  curvature  of  the 
detonation  front. 

2.  Before  proceeding  with  a  general  analysis  of  the  curved- 
front  theory,  let  us  derive  the  detonation  velocity  of  an  ideal  wave 
in  an  imperfect  gas,  assuming  an  equation  of  state  of  the  form 

pv  -  nRTex  (2.1) 

where 

x  -  K/v.  (2.2) 

Previous  studies  have  used  a  constant  co-vo lurae  equation  of  state. 

This  derivation  results  from  a  simultaneous  solution  of  the 
hydrodynamic  equations.  The  equation  of  continuity  for  a  plane 
wave  is 

U/D  -  v/vo  (2.3) 

and  the  equation  of  motion  is 

p  *  U2/v  »  D2/vo  +  p0  (2.4) 

Furthermore,  the  equation  of  energy  conservation  is 

AE  ♦  pv  ♦  1/2  U2  ■  4Q  ♦  1/2  D2  ♦  p^ 


(2.5) 


where 
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AE  =  Cy  (T  -  Tq)  =  Cv  (pv/nReX  -  Tq) 


(2.6) 


j't  cvdT 


and  Cv  =  xTTY"  9  since  the  internal  pressure  p^  «  T  -p  is 
o  >  /v 


zero  by  equation  (2,1). 
Substituting  (2.6)  into  (2.5), 


where 


L'pv  +  1/2  (U2-D2)  -  Q 


Q  »  4Q  +  CyTo 


(2.7) 


(2.8a) 


L'  -  C^tiRe*  +  1 

and  where  pQ  is  assumed  to  be  negligible  compared  to  p.  The 
Chapman  Jouguet  condition  is 

°  ■  • 

Differentiating  (2.1)  with  respect  to  v,  and  using  the  thermo¬ 
dynamic  relations 


(2.8b) 


(2.9) 


we  have 


where 


Hence,  by  (2.9), 


(SL)  .  _  E_  (to)  ..x 

W; S  Cv  Wj  v 

.  _  £  1 

Vdv/s  V  L 

1/1  ■  1  ♦  nRe^C  ♦  x 


D2/vo2  -  p/rL 


(2.10) 


(2.11) 


Fro*  (2.3),  (2.4)  and  (2.11) 


D2/p  *  vo/(l-v/vo)  -  vo2/vL 
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or 

v/vq  =  l/(l+L)  (2.12) 

For  the  explosive  EDNA,  the  ratio  of  v/vq  given  by  (2.12) 
has  been  calculated  as  follows  for  x  in  the  range  0.8  to  2.2, 
using  Cy  and  n  data  obtained  in  Technical  Report  No.  VI. 

x  0.8  1.0  1.2  1.6  1.6  2.8  2.0  2.2 

v/vq  0.708  0  .  724  0.745  0.763  0.778  0.792  0.805  0.815 

This  illustration  indicates  the  degree  of  error  introduced  by  the 
approximation 

v/vq  «  3/4  ,  .(2.13) 

which  will  be  assumed  reasonable  for  most  solid  explosives 
and  used  in  this  discussion. 

Also,  from  2.7,  using  2.3,  2.11  and  2.12, 

(L'L  ♦  1/2)  (1+L)"2  -  1/2  •  Q/Dj2 

where  is  the  ideal  detonation  velocity. 

3.  Let  us  know  proceed  to  analyze  the  non-ideal  detonation 
wave  under  the  assumptions  of  the  curved-front  theory.  Assume  as 
a  first-order  approximation  that  any  small  region  of  the  wave-front 
is  spherical.  Then  the  normal  detonation  velocity  D  may  be  cal¬ 
culated  from  the  conservation  equations  appropriate  to  a  spherical 

detonation  wave  with  r  the  local  radius  of  curvature. 

o 

If  a  is  the  angle  the  normal  makes  with  the  axis  of  propagation 
and  if  the  wave  proceeds  without  a  change  of  shape  (steady-state) 

D-D  cos  a 
o 

where  Dq  is  the  propagation  velocity  of  the  wave  itself. 

Integration  of  the  differential  equation  of  continuity  for  a 
spherical  shock  leads  to 

U/D  -  vS/vo 


(3.1) 


-  5  - 
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By  (3.5),  (3.6)  and  (3.8), 


2L+(i+  e2) 


If  Now  we  set 


2<p  »  1  -  9  2, 
then  (3.9)  becomes 

b  *  ihr  • 

By  (3.8)  and  (3.11), 

z  -  - j 

(1+L^)2 


Also, 


2  l-2(p 

w  ■  - x—0 

(1+1-9) 

Substituting  (3.12)  and  (3.13)  into  (3.7), 


(3.9) 

(3.10) 

(3.11) 

(3.12) 

(3.13) 


(U«  ♦  1/2  (1-2*))  - i— K  -  1/2  -  %2  , 

(1+L^r 

or 

A  ♦  L»  -  1/2  92  -  (1+1-9)2  %2  ,  (3.14) 

where 

A  •  L  (L'  -  1  -  1/2  L). 

Under  ideal  conditions,  9  -  0,  i.e., 

A  -  (1*L)2  Q/Dt2 


(3.15) 
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Hence,  from  (3.14), 


1 


iilkli2  Ok  , 

(l+L)2  D2 


(3.16) 


For  small  values  of  cp,  a  good  approximation  to  (3.16)  is 


D_  l-M'g 
"  1+N«9 


(3.17) 


where 

M'  -  l/(l+L),  N'  -  2(L^y.L  (3.18) 

As  an  example  of  the  percentage  error  involved  in  the  approxi¬ 
mation  of  equation  (3.17),  consider  the  case  of  x  *  1.8,  Cv  •  0.44, 

Cv  *  0.38  and  n  ■  34.94  for  the  composition  EDNA.  The  error  is  less 

than  2  per  cent  for  <p  less  than  0.1,  which  corresponds  to 
greater  than  0.9  (and  less,  of  course,  than  1.0). 

4.  In  order  to  determine  the  shape  of  the  detonation  front, 
it  is  sufficient  to  find  a  parametric  representation  of  the  cross- 
section  of  the  front,  recognizing  the  symmetry  of  the  wave  around 
the  axis  of  the  charge  (see  Fig.  1).  For  any  point  of  the  front,  the 
radius  of  curvature  is  given  by  rQ  ■  ds/da,  where  s  is  the  distance 
measured  along  the  front,  and  a  the  angle  the  normal  makes  with  the 
wave-front.  Also,  from  (3.3)  and  (3.10), 


*  -  -  2/3  in  (l  -^)  (l  ♦  1/3  In  (l  -  |j|J 

A  plot  of  this  function  (rig.  2),  gives  a  good  linear  approxima¬ 
tion  in  the  form 

*  •  0.7  a Jr  (4.1) 

o  o 


for  a  /r  <  0.75 
o  o 


Equations  (4.3)  give  the  parametric  representation  of  the  wave 
front.  For  the  case  of  a  solid  explosive  surrounded  by  air,  the 
final  angle  can  be  taken  to  the  tt/2.  For  this  case  we  have, 


yc/a  =  (M+N)Kc  -  N  (4.5) 

where 


The  universal  function  K  is  given  in  the  following  table. 

C 

6  K  6  K 

c  c 


0.000 

1.00 

.850 

3.93 

.200 

1.19 

v  875 

4.44 

.300 

1.32 

.900 

5.11 

.400 

1.45 

.925 

6.13 

.500 

1.69 

.950 

7.66 

.600 

1.99 

.975 

11.86 

.700 

2.45 

.990 

19.90 

.750 

2.78 

*994 

26.32 

.000 

3.24 

.997 

3?.  12 

-  9  - 


As  an  illustration,  using  theoretical  data  for  the  composi¬ 
tion  EDNA  as  given  in  Technical  Report  No.  VI,  the  following 
results  are  obtained.  Figure  3  is  a  plot  of  the  ratios  Dq/D^ 


& 


to  a  /y  for  various  densities, 
o  c 


1.58  a  /r 
1.20  0  0 
0.93 


0.997  0.994  0.990  0.975  0.950  0.900  0.850  0.750 

0.022  0.032  0.042  0.072  0.11  0.18  0.24  0.36 

0.025  0.037  0.049  .084  0.13  .21  0.28  0.42 

0.028  0.041  0.055  .094  0.15  .23  0.31  O.46 


The  relationship  of  D  /D.  to  a  /y  is  often  approximated 

0  1  o  c 


linearly.  Such  an  approximation  introduces  about  a  five  per  cent 
error  for  ideal  conditions,  as  can  be  seen  by  a  linear  continuation 
of  the  curves  in  Figure  3. 


